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Phenomenological scaling arguments suggest the existence of universal amplitudes in the finite- 
size scaling of certain correlation lengths in strongly anisotropic or dynamical phase transitions. For 
equilibrium systems, provided that translation invariance and hyperscaling are valid, the Privman- 
Fisher scaling form of isotropic equilibrium phase transitions is readily generalized. For non- 
equilibrium systems, universality is shown analytically for directed percolation and is tested nu- 
merically in the annihilation-coagulation model and in the pair contact process with diffusion. In 
these models, for both periodic and free boundary conditions, the universality of the finite-size 
scaling amplitude of the leading relaxation time is checked. Amplitude universality reveals strong 
transient effects along the active-inactive transition line in the pair contact process. 



The notions of scaling and universality are central to the modern understanding of critical phenomena, see e.g. 
■ Besides the well-known universality of the critical exponents, universality is also manifest for many critical 
amplitudes, as reviewed in B. Here we are concerned with universal amplitudes which arise in finite-size scaling. For 
example, consider a statistical system at equilibrium, e.g. a simple ferromagnet, which is described by an isotropic 
and translation- invariant field theory in its continuum limit. Close to its critical point, and on a lattice with finite 
extent L, Privman and Fisher [^| showed that in any dimension d below the upper critical dimension d* the singular 
part of the free energy density / and the inverse correlation lengths £j satisfy the scaling form 



where z\ — CitL 1 ^ , z 2 = CihL^^ 1 ' '< lv are the scaling variables, the reduced temperature t = (T — T c )/T c and the 
reduced magnetic field h — H/T c describe the distance from the critical point and /?, 7, v are the standard equilibrium 
critical exponents. The index i distinguishes different correlation lengths, e.g. i = a for the spin-spin correlation 
length or i — e for the energy-energy correlation length in a simple ferromagnet. Furthermore, the scaling functions 
Y and Si are universal functions and all non-universal properties of a given model are condensed in the values of the 
non-universal metric factors C\ and C 2 - Although the functions Y, Si are universal, they do depend on the boundary 
conditions and of the geometric shape of the lattices under study. 

It follows from (|l|) that if the model parameters are tuned to t = h = so that the model is at its bulk critical point, 
the finite-size scaling of free energy density / = L~ d Y(0, 0) and of the inverse correlation lengths i^ 1 — L~ 1 Si(Q, 0) 
is described in terms of universal finite-size scaling amplitudes. This qualitative statement can be made quantitative 
in 2D for systems defined on an infinitely long strip of finite width L through a by now classic conformal invariance 
argument which for example explicitly relates the 5^(0,0) to known critical exponents |(|. Concrete model studies 
have confirmed this many times, as reviewed in More recently, similar relations have been conjectured also in 

3-D from the results of numerical studies in toruslike |7]|| and spherical geometries |]|l(| ■ 

When considering the situation of dynamical scaling, where time and space scale differently, the extra degree of 
freedom might appear to exclude the presence of universal finite-size scaling ampltitudes in the sense of eq. (|l|). Here, 
we shall ask under what conditions the arguments of Privman and Fisher |Bj can be generalised also to anisotropic 
scaling. Specifically, for a system undergoing a phase transition with anisotropic scaling which is defined on a lattice 
with finite extent L in the spatial direction but is infinite in the temporal direction. Then (under conditions to be 
detailed in sections 2 and 3) the spatial correlation lengths should satisfy the scaling form 
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I. INTRODUCTION 



f{t,h)=L- d Y{ Zll z 2 ) , C rl =L- 1 S i (z 1 ,z 2 ), 



(1) 




(2) 



in a notation analogous to (|l|) and where again the Si are universal functions. 
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While the generalization to equilibrium anisotropic scaling is a rather straightforward extension of the methods 
valid for the isotropic case S, the non-equilibrium situation is more difficult. We shall show how correlation length 
amplitude universality can be established for systems in the directed percolation universality class. Furthermore, 
numerical data from some reaction-diffusion models are also in agreement with this and suggest that the scaling form 
(||) might be generally valid. 

This paper is organised as follows: in section 2, we present the scaling arguments leading to the recognition 
of universal amplitudes for anisotropic scaling in equilibrium. In section 3 non-equilibrium dynamical scaling is 
discussed, with emphasis on the directed percolation universality class. In section 4, we test the amplitude universality 
by analysing the finite-size scaling of the leading inverse relaxation time in the annihilation-coagulation model and 
the pair contact process with single particle diffusion. In section 5 we give our conclusions. An appendix contains the 
Bethe ansatz calculation of the relaxation time in an analytically solvable special case. 



II. UNIVERSAL FINITE-SIZE AMPLITUDES IN EQUILIBRIUM 

Consider an equilibrium system with anisotropic scaling in two distinct directions. Systems of this kind arc known 
since a long time, e.g. in Lifshitz points Jll]] or anisotropic uniaxial magnets see [jl3| p~6[| for reviews or else in 
quantum phase transitions jf?]], see (T^| for a recent book. One can derive the phenomcnological scaling of the physical 
observables from the covariance of the correlators under scale transformations and reconstruct the thermodynamics 
this way, see ||,|| . In line with conventional phenomenological treatments of strongly anisotropic scaling [|l9|,^0| , we 
assume for the two-point functions the scaling form 

Gi (r x , ry ;t,h)= b~ 2x < G< , £J; tfi* , ) (3) 

where b is the rescaling factor, t, h refer to physical quantities like the reduced temperature and the reduced magnetic 
field, 9 is the anisotropy exponent and yt,yn, Xi are scaling exponents. The index i refers to different physical quantities 
of which the two-point function is formed, e.g. i = a for the spin operator or i = e for the energy density (for simplicity, 
we use throughout a notation analogous to simple ferromagnets) . At criticality, t = h = 0, and one has 

GUrx.O) ~ r-I 2 ^ , G ff (0,r||)~r- 2aW9 , 

G £ (r ± ,0)~r^ ; GMr\\)~r- 2x ' /e . (4) 

For a strongly anisotropic equilibrium system, r± and rii correspond to different directions in space. This case is for 
example realised at the Lifshitz point in spin systems with competing interactions like the ANNNI model |l5|,|l| • For 
brevity, we shall refer to the directions r± as 'spatial' and to the directions ni as 'temporal'. 
We use the equilibrium fluctuation-dissipation theorem 

X = fdr ll d d r ± G a (r ± , n ) , C= f dr\\d d r± G £ (r±, r\\), (5) 



where d is the number of 'spatial' dimensions and x, C are the susceptibility and specific heat. We shall work with a 
single 'temporal' direction throughout, but generalisations are obvious. Units are such that the critical temperature 
T c = 1. From (||) and integrating, one gets immediately the scaling forms for \ = x(^^) = —d 2 f/dh 2 and for 
C = C(t,h) = —d 2 f/dt 2 . Here / is the (singular part) free energy density. From the scaling of x an d C, it should 
satisfy the scaling 

f(t,h) = b 0+d - 2x "- 2yh f{tb y \hb Vh ) = b e+d ~ 2x ^ 2vt f{tb y \hb Vh ). (6) 

These two forms are consistent, if x a + yh = x e + yt — w. In fact, the above argument works for any pair of 
scaling operators like a, £ and their conjugated scaling fields h, t. Therefore, the value w must be independent of 
all physical scaling operators which might be present in a given model. Next, we define the standard static critical 
exponents a, f3, 7 as usual, e.g. Q. Also, out of criticality, we expect an exponential decrease of the two-point function, 
characterised by the correlation lengths £j_ n ~ i _!/J - ii (at h = 0). From this, we find that u± = l/yt, 8 = v\\/v±_ and 
f3 + 7 = yh/yt- From dimensional counting, we expect w = d + 9 and the free energy density then scales as 

f(t,h) = b- d - e f(tby\hby-). (7) 
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As usual, we explicitly assume the absence of dangerously irrelevant scaling fields EIMO]. We then recover the 
hyper scaling relation 2 — a = v\\ + di/±. The common scaling form for the thermodynamics is found by scaling out b 
and introducing the conventional critical exponents: 

f(t,h) = A 1 \t\ 2 - a W ± (A 2 h\t\-^), (8) 

where W ± are the universal scaling functions which are obtained for t > and t < 0, respectively, and A\ i2 are 
non-universal metric constants. At this level, the anisotropy of the scaling of the two-point function only appears in 
the generalized form of the hyperscaling relation. 

Equation (g) will be the starting point for our discussion of finite-size scaling. In what follows, we consider a 
situation where the 'spatial' directions are of finite extent L whereas the 'temporal' direction remains infinite. In the 
same spirit as Privman and Fisher [||, we assume that the finite-size scaling behaviour is governed by the 'spatial' 
correlation length £j_ only and we write 

f(t, h; L) = A 1 \t\ 2 - a W ± (A 2 h\t\-^; L^ 1 ) , (9) 

where the bulk 'spatial' correlation length £j_ = £ot~ l/± - Note that there is no extra metric factor in the second 
argument of W^, whereas £o is non- universal. To simplify the notation, we assume that there is no phase transition 
for L finite, but this restriction could be removed analogously to the equilibrium case |^| . We emphasize that the 
'temporal' direction remains infinite, otherwise we would have to deal with two distinct finite length scales. 

Following the ideas developed by Privman and Fisher we have to trace the non-universal constants, taking into 
account the anisotropic scaling. For that, it is sufficient to study the 'spatially' infinite system. We expect the scaling 
form of the connected spin-spin correlator, see also (|), 

G a (r±,r r ,t,h) =D D 1 r- 2x "X ± {r^/^,D () r ] \/^-D 2 h\t\-P-<), (10) 

where A^ is a universal scaling function and -Do, 1,2 are non- universal metric factors. From the fluctuation-dissipation 
theorem (||) one has 

X (t, h) = D i e ± /v± X ± (D 2 h\t\-^) . (11) 

where is a new scaling function obtained from X^. Now, for the non-connected spin-spin correlator, one has in 
the same way, introducing new universal scaling functions 

T ff (r x ,r\\;t, h) = D Q D 1 rl 2x 'Z ± (r ± /£ ± , D r\\/&; D 2 h\t\-^) (12) 

At this point, we assume translation invariance with respect to r± and rn. Therefore, there should exist a mean 
magnetization m which is independent of r±_ and rn and which can be found considering r CT at large separations 
r±,rf 

r„(rx,r\\;t,h) = DqD^'Z^ (l, A)r-||A±; D 2 h\t\~^) (13) 
I i 

m 2 (t, h) = B Q D 1 C 2x °Z ± (D 2 h\t\-^) , (14) 



where the arrow indicates taking the limit of large 'spatio-temporal' separations. Because of translation invariance, 
m 2 should become independent of Dq. On the other hand, applying standard thermodynamics to the free energy (|^) 
yields 

m(t,h)=A 1 A 2 \tfW± (A 2 h\t\-P-~>) , (15) 
X (t, h) =A 1 A 2 2 \tr\V 2 ± (A 2 h\t\-P-"<) , (16) 

where W±(x) = d n W ± (x)/dx n . 



xVd 
& 



We now compare eqs. ( |1 1| ) and (|16[). Letting first h = 0, we find 

=A 1 A 2 U 1 . (17) 
Comparing the arguments of the scaling functions, we have 
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D 2 = A 2 U 2 . 



(18) 



Next, we compare eqs. (|IJ) and (|l5|) and find for h = that 

D D^~ 2 ^ = AfA 2 2 U 3 (19) 

(since x a = fl/v±). Here, f/1,2,3 are universal constants whose universality follows from the universality of the scaling 
functions considered. Using the hyperscaling relation 7 + 2(3 = (d + 8)v±, we find that 

Ai^ +e D^ = Q x = Ux/Us 

D2A2 1 =Q 2 = U 2 (20) 
D iJ(u±(d+e)) D ^ A -i- 1 /(u ± (d+e)) A - 2 _ g 3 _ u i- 1 /(v ± (d+e)) u ~t/(v ± (d+e)) 

and the Qi,2,3 are universal constants. 

Finally, we come back to the finite-size scaling behaviour. In eq. (^), we replace £0 by A\ using (pc|). Scaling out 
L and using again hyperscaling, it is easy to arrive at the scaling form 

f(t,h;L) = L- d - D o Y (c l tL 1 /^,C 2 hL {f}+ ' l) l w ^ , (21) 

where Y is a universal scaling function and C\ i2 are non-universal metric factors related to A\. 2 . In contrast with 
the isotropic situation ([l]), we see that the finite-size scaling amplitude of the free energy is no longer universal. 
Furthermore, since £|| = ^/Z?o, we have 

/(i,0;£)£i(t,0;Z)£||(t,0;L) = D^f{t, 0; L)^[ +9 (t, 0; L) -> univ. constant, (22) 

which holds because of (^0[). Therefore, we expect for the 'spatial' correlation length 

^J 1 = L~ 1 S (c x tL l / Vl - , C 2 hL^ + ^ u A (23) 

with a universal scaling function 5* and the same metric factors C\^ 2 as in (pll). While this analysis was phrased 
in terms of the transverse spin-spin correlation length £± = £j_,<7-j similar arguments should hold for the 'spatial' 
correlation lengths £x,j °f an y other physical observable, with S in (|2^) being replaced by an appropriate function Si. 
The scaling functions Y and Si should depend on the boundary conditions and, if d > 2, on the shape of the finite 
'spatial' domain. Note that L refers here to the physical length, which cannot be equated to the number of sites N 
times the lattice constant in non-square lattices pP^|. For a recent example of this in the 2D Ising model context, 
see @. 

Finally, the 'temporal' correlation lengths 6u should read 

^ = L~ e D R t (CttL 1 /^ , CthL^y"^ (24) 

with universal scaling functions Ri and again the same metric factors C\ i2 as before. The value of Dq is related to the 
dimensionful constant which occurs in the energy-momentum dispersion relation of the underlying continuum field 
theory and cannot be found straightforwardly. However, at criticality (t = h = 0) ratios of 'temporal' correlation 
lengths ClLi/Cllj should tend to universal constants in the L — > 00 limit. The universality of these ratios would not 
have immediately been obvious from straightforward finite-size scaling. 

Eq. ( p3|) is the main result of this section. It provides the natural generalisation of the Privman-Fisher form (|l|) to 
the case of equilibrium anisotropic scaling. It is immediate to include further physical parameters into the analysis. 
We emphasize that (i) translation invariance and (ii) hyperscaling was required in deriving this result. We stress that 
we considered finite sizes in the 'spatial' direction and obtain universality for the spatial correlation lengths £1^ only. 



III. UNIVERSAL FINITE-SIZE AMPLITUDES OUT OF EQUILIBRIUM 

We have seen that in equilibrium systems with anisotropic scaling and in a geometry where the 'spatial' directions 
are of finite extent L while the 'temporal' directionis infinite, the 'spatial' correlation lengths have a universal finite- 
size scaling amplitude. We now ask whether this result generalizes towards more general forms of dynamical scaling, 
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without appealing to the special properties of equilibrium systems. Fluctuations in non-equilibrium systems can be 
treated in terms of dynamic functionals via Martin-Siggia-Rose theory p^-|2?[. To be specific, we shall work in a 
setting of reaction-diffusion processes, of which directed percolation is a common example, see [^8] [HJ for reviews. 
We shall continue to denote time by rn and space by r±. As before, t measures the distance from the steady-state 
critical point and h denotes an external field (e.g., for directed percolation t = p — p c and h is the rate of a process 
— > A). For the sake of technical simplicity, we shall assume translation invariance throughout. As in equilibrium, 
we have to trace the non-universal metric factors and this is most conveniently done in the bulk. 

Physical quantities of interest are the mean particle density p, the survival probability P and the pair connectedness 
function G = G(r'j_,r'prx>r\\), which is defined as the probability that the sites (?"j_> r n) anc ^ ( r -L' r ||) are connected 
by a direct path |2j],|3(J. Because of translation invariance G = G(r' ± — r±, ri — r\\), which will be used throughout. 
These quantities are expected to satisfy the scaling behaviour 

P (r ± ,r f t } h) = b-^p r -l;t h y^hby-) = d 1p ^ x "s ± A>^; D 2 h\t\-^ 

P(r x ,r r ,t, h) = b~ xp P jljt&w hb^) = D 1P (^±,D ^;D 2 h\t\-v^ (25) 

G(r±, ry ;t,h) = 6-°G (^, ^; W , hb^) = D 1G C i XG Q ± (j^, D ^ D 2 h\t\-y»'*) 

where the x's are scaling dimensions and y^h renormalization group eigenvalues, the D's are non-universal metric 
factors, £, T, Q are universal scaling functions where the index distinguishes between the cases t > and t < 0, 
£j_ = Co |*| ~ v± is the spatial, £n = ^/Dq is the temporal correlation length and z is the dynamical exponent (as before 

y t = l/v±). 

In the steady state, and for h = 0, one expects p ~ t 13 and P ~ . In general, the two exponents f3 and (3' are 
distinct from each other. For spatial translation invariance, the dependence on r± drops out for both p and P and in 
the steady state (i.e. rn — > oo) one has 

P (t,h) = Dxp^'^ (p 2 h\t\-™/*) \tf 

P(t, h) = D 1P ^' , ^F ± (D 2 h\t\-y^) \tf (26) 



where x p = 0/u±, xp = /3'/i>± and 8 ± — lim r ,|_,oo £^ and similarly for T . We also consider the auto-connectedness 
(that is r± = rj_) in the steady state 

G(0, oo; t, h) =: G(t, h) = D 1P ^*°Q± (p 2 h\t\-**l*) \t\ xa "± (27) 

In the active phase (t > 0), the surviving clusters will create an average density ~ \tf in the interior of the spreading 
cone. Therefore, the auto-connectedness should in the steady state saturate at the value 32 

G(t,h) = p(t,h)P(t,h) (28) 

Comparison of the scaling forms then yields, setting h — 0, 

x G = ifl + P)fvx , Diq = Di p Dip (29 ) 

y (0) 

Usually, xq = d — 6z is expressed in terms of the initial critical slip exponent 6 [[33|, which makes it apparent that 
the expression (p9| ) is in fact a generalized hyperscaling relation |2^,^,^2| . 
Next, we consider the total mass M of the cluster [f30|, given by 

M(t,h) := [ d d r ± /""dry G(r ± ,r f t 7 h) = ^-^/"^ ^D 2 h\t\- y ^ y ^ (30) 



Ei' 1 



T=T± 



where eq. (25) was used and is a new universal function related to Q ± . Also 

7 = du x + v\\ - & - & (31) 
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which is the analogue of the hyperscaling relation of the equilibrium systems. 

While the discussion so far has been completely general, we now appeal to two properties which are valid for 
systems in the directed percolation universality class, but need not be generically valid. First, we consider a directed 
percolation process in the presence of a weak field h (physically, h parametrises the rate of a particle creation process 
— ► A). A site at a given time becomes active if it is connected with at least one active site in the past, where a 
particle was created by the field. The number of such sites is equal to the cluster size, the probability to become 
active is given by the density J3(| 

p(t, h) ~ 1 - (1 - h) M{t ^ ~ hM(t, h) (32) 
dp(t, h) 



for h small. Therefore, 



M(t,0) 



(33) 

h=Q 



dh 

Comparison with the scaling forms for p and M leads to 

Vh/Vt = P + 7 , D 1P = D D 2 ^ P+ ^A ± (34) 

where is an universal amplitude. Second, directed percolation is special in the sense that there is a 'duality' 
symmetry which can be used to show that |p4| 

p(t,h) = P(t,h) (35) 



As a consequence, /? = /?' and D\ p = Dip for directed percolation and we thus have, combining eqs. (|26|,[30 34) 

p(t,h) = D Q D 2 ^ d - z \tfMf (D 2 h\t\-^) 
M(t, h) = D D%£o d - z \t\-"<Mt (D 2 h\t\-P-~<) (36) 



with universal functions AA^(x) = d n JA ± (x)/dx n and where the hyperscaling relation eq. ( pl| ) has been used. We 
therefore recover the analogues of eqs. ([l5 16) found in equilibrium. Finally, we define a new function p = p{t,h) by 
p(t, h) = dfi(t, h)/dh, which implies 

/i(t, h) = D ^ d - Z \t\ {d+Z ^M ± (D^-P- 1 ) (37) 

In particular, as we did before at equilibrium, we have because of £|| = ^j_/Dq that 

H(t, 0)£l (t, 0)^|| (t, 0) -» iiiiiv. constant, (38) 

which is indeed the analogue of the result ( f22| ) . 

At last, we consider a geometry of finite size L in space but of infinite extent in time. Again, by analogy with 
section 2, we postulate that in this finite geometry merely the scaling functions are modified 

Mt = Mi{D 2 h\t\-^;Lil 1 ) (39) 



and without introducing any further metric factor. Indeed, we can then scale out L and, because of eq. (|38[), arrive at 
the same scaling forms ( p3| , p4[ ) as had been found before for the spatial and temporal correlation lengths in anisotropic 
equilibrium systems, at least for systems in the directed percolation universality class. Since directed percolation is 
known to be equivalent to an equilibrium (in fact, purely geometrical) problem, the existence of universal finite-size 
amplitudes in this class is not too surprising and might have been anticipated from the discussion in the previous 
section. 

While the spatial correlation length £j_ may not be always a very accessible quantity, its universality may also be 
tested by considering the spatial moment (let t = h = for simplicity) 



R { l> ■■= (rl) 



SA(L)^r±S"<ir\\rlG(r ± ,r r ,L/a) 
W^JTdl G(r x ,r,|;L/a) 
W)dV gdni rr' g g ± (rx,r||;L/q) 
/acl/o^x /~dr|| rI-"0±(r X) r„ ; L/U) 
GH„(i/ex) 

L n E n (L/U) (40) 
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where A(L) is a d-dimensional hypercube of linear extent L and S„ and S„ are universal functions. Since there is 
no metric factor in the argument of S n , the universality of the finite-size scaling amplitude of £j_ is equivalent to the 
universality of the finite-size scaling amplitude of Ii£ (on the other hand, the temporal moment (r») ~ {D^L 2 ) 11 has 
a non- universal amplitude). The universality of these moments is a somewhat stronger statement than the universality 
of certain ratios of moments (p n ) of the particle density p which has recently been verified in ID and in 2D for several 
models in the directed percolation universality class ]35| ]. 

In summary, we have seen that for systems in the directed percolation universality class, the special properties 
eqs. ( |33|,35| ), taken together with the general relation (|2|), are sufficent to rederive the universal finite-size scaling 
form ^22) of the spatial correlation length, in spite of the absence of the fluctuation-dissipation relation. It is not yet 
clear wether there exist more general arguments which would permit us to arrive at the same result without appealing 
to either (|3^) or (|35|). However, we shall in the next section present numerical evidence that the universal finite-size 
scaling forms (|2^) for £j_ or ( pi| ) for £|| might be more generally valid. 



IV. REACTION-DIFFUSION PROCESSES 

The new information contained in (|23| ) which goes beyond the standard RG ideas is the universality of the finite-size 
scaling amplitude Lt;^ 1 precisely at criticality (and similarly the universality of all ratios £||,-t/£|| ,,?)• Since for the time 
being, the derivation of this universality for non-equilibrium systems appears to be restricted to directed percolation, 
we use the pair contact process and the annihilation-coagulation model to test the universality hypothesis advanced 
in sections 2 and 3 quantitatively. 

The pair contact process []36f has been intensively studied recently. It is a reaction-diffusion system, where particles 
move and react on a lattice. While each lattice site can be either empty or be occupied by a single particle, the 
following microscopic moves are permitted 
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rate 


d 





(41) 

and are parametrized by the diffusion constant d and the pair annihilation rate p. 

While in the case without diffusion (d = 0), the steady-state transition between the active and the absorbing state 
was found to be in the directed percolation universality class |36|]35|^^| , the effects of adding diffusion were first 
studied using field-theoretical methods, considering a bosonic field theory without any restriction on the number of 
particles per site (which leads to a divergent particle density in the active phase) |^8[ . It was shown that the entire 
absorbing phase is critical and in the universality class of diffusion- annihilation (see below). Because of the non- 
renormalizability of the underlying field theory, no quantitative information about the transition towards the active 
state could be obtained. The first quantitative informations were obtained |39| through the use of density matrix 



renormalization group (dmrg) techniques [}40|-|42| . The steady state phase diagram is shown in figure [T| an d there is 
a general agreement between dmrg and Monte Carlo studies on the location of the critical point p c (d) |39| , ^3 -46|. 
The annihilation-coagulation model is formulated in the same way, with the allowed reactions 

AA -> A$ , $A with rate dj 

AA -» 00 with rate Ida (42) 

A® <-> 0A with rate d 

parametrized by a and 7 for annihilation and coagulation, respectively. The long-time behaviour is the model is 



always algebraic, i.e. the mean particle density p(t) ~ t 1//2 in ID, sec [£8-31 47| and references therein 
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FIG. 1. Steady state phase diagram of the pair contact process. The dotted lines are the phase boundaries according to pair 
mean field theory, while the full line gives the active-inactive transition in ID. DP marks the steady-state transition in the 
directed percolation universality class at d = 0. 



In this paper, we shall use the Hamiltonian formulation of reaction-diffusion processes 48-5^], which starts from 
the master equation 



d\P(t)) 
dt 



-H\P(t)) 



where \P(t)} a state vector and H is referred to as "quantum" Hamiltonian (for recent reviews, see 
chain with L sites, H is a stochastic 2 L x 2 L matrix with elements 



(a\H\r) 



-w 



(43) 
For a 

(44) 



where |er), |r) are the state vectors of the particle configurations <t, r and w are the transition rates. It is well-known 
that the ground state energy of the pair contact process E = Ei = is twofold degenerate for ^ |39[] . The 
energy gap r = E2 — Eg, calculated from the second excited state, is the inverse relaxation time or, in the notation 
of section 2, the inverse 'temporal' correlation length £,7 = T towards the steady state. We shall consider both free 
and periodic boundary conditions. 

1. First, we discuss the pair contact process (ffl|). We find the following, surprisingly simple, finite-size scaling 
behaviour for the gap Tl in the entire absorbing phase, that is, for all p > p c {d), namely 



where a depends on the boundary condition 



2ir 2 

_2 



; periodic 
; free 



(45) 



(46) 



but is independent of both p and d. 

Before deriving (46), we argue that this result confirms the universality of the correlation length amplitudes dis- 
cussed in section 3. Indeed, for systems in the diffusion-annihilation universality class already the inverse 'temporal' 
correlation length T has a universal amplitude, provided that the value of the diffusion constant d is fixed. The 
universality of the reaction-diffusion process 2 A — > has been discussed using field theory methods 53 - 55 1 . From 
the renormalisation group equations, it can be shown that the value of the diffusion constant d is not renormalised 
through the effects of the interaction of the particles and simply stays at its bare value. The bare value of d is the 
value it has in the original lattice formulation of the problem. Since the diffusion constant sets the time scale, we 
expect for the gap T ~ L~ 2 d. These calculations |53|-|5f| apply to the process 2A — > which corresponds to the 
case p = 1 in the model at hand. However, it is known that in the entire absorbing phase, the extra interactions 
coming from the reaction 2 A — > 3A are irrelevant Therefore, they should not modify the value of d (since we 
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only consider here the inactive phase, we leave aside the question how d evolves under renormalization in the active 
phase). Consequently, the proportionality of T and d in eq. ( [45| ) comes from the non-renormalisation of d. That 
non-renormalisation is a special property of the diffusion-annihilation universality class. That is consistent with the 
scaling form ( p^ ) for £|| and we can identify Dq — d. Finally, the p-independence of the amplitude a in eq. ( |46| ) 
is an example of the universality of the finite-size scaling amplitude R(0, 0) in eq. (p4|), which in turn implies the 
universality of S(0,0) in eq. (p3|). 

We now derive eqs. ( p^ , [46| ). We need the lowest non- vanishing eigenvalue r = E2 of the quantum Hamiltonian H. 
For p = 1, it turns out that the spectrum of H is equal to the spectrum of an XXZ Heisenberg quantum chain -ffxxz- 
The lowest gap of i/xxz can be found from the coordinate Bethe ansatz (5^J|^] . This reproduces eq. ( f46| ) for p = 1 
and all values of d. The details of the calculation are presented in the appendix. 

We point out, however, that for d = 1/2 and free boundary conditions only, the Bethe ansatz equations have a 
closed-form solution. The exact lowest gap for any finite number of sites L is 

7T 1 

T L = 1- cos j— ; d = - , p = 1 , free b.c. (47) 

in agreement with ( |45| , ^6| ). This had been conjectured before on the basis of numerical data (closed- form solutions 
for slightly different A = 1/2 XXZ chains have been discussed recently in p8[). 

For p ^ 1, the PCPD is not related to any known integrable model and we revert to numerical methods. We 
consider the normalized amplitudes 

A ( P = L 2 T L /(2n 2 d) , A ( P = L 2 T L /(n 2 d) (48) 

defined for periodic (P) and free (F) boundary conditions, respectively. If and only if ( p6[ ) is correct, the amplitudes 
A [ P and A^ ] should converge towards unity in the L — > 00 limit. 



In table |I| we show data for periodic boundary conditions. These were obtained from diagonalizing H through 
the standard Arnoldi algorithm Q. Translation and parity invariance were used to blockdiagonalize H, and lead to 
matrices of size ~ 2 L / (2L) for L sites. In table |n| data for free boundary conditions are shown. They were obtained 
by applying the density matrix renormalization group (dmrg) method to the pair contact process [j39| . 

(P) (F) 

Clearly, the data for both A y L ' and A L ' at finite values of L are in general quite far away from unity. We also 
see that the raw data tend to be closer to unity for smaller values of d. However, since the systematic variation 
of these amplitudes with L is huge, a precise L — > 00 extrapolation must be performed. We have used the BST 
extrapolation algorithm [^9| which has established itself as a reliable and precise method for the extrapolation of 
finite-lattice sequences arising in both equilibrium and non-equilibrium critical phenomena. The parameter oj describes 
the (effective) leading finite-size correction of a sequence Al = Aoo + AiL^^ + . . . and must be chosen to optimize 
convergence |j|j59|. 

In all cases, we find that the extrapolated amplitudes 

Ag> := lim 4 P) ~ 1 , Ag> := lim A<?> ~ 1 (49) 

L — *oo L — >oo 

within the numerical accuracy of the extrapolation and in full agreement with eq. ( f46|) . The need for L — > 00 extrap- 
olation also means that the universal infinite-size amplitudes may be hard to see in, say, Monte Carlo simulations. 
We also give the effective value of u> for each sequence. In view of the exact result ui = 1 for p = 1, see (|5^^), it is 
satisfying that u stays close to one. 

In general, for a given d convergence is best for p close to 1 and decreases when p is lowered. The lowest values of 
p given in table [| are more or less the smallest ones for which a reliable convergence of the amplitudes could still be 
observed. In varying d, we see that data converge best for relatively small d and that close to d = 1 the cross-over 
towards mean-field behaviour [p9[ affects the finite-size scaling of the amplitudes. Comparing the data for periodic 

— (F) 

and free boundary conditions for the same values of p and d, we observe that the A^ are closer to unity than the 
A)x,' for the same value of L (and in contrast to the usual expectation that finite-size corrections should be smallest 
for periodic boundary conditions). For p = 1, these remarks are confirmed analytically, see eqs. j57|j6^ ). 

All in all, the extrapolated amplitudes converge over a large range of values of p and d towards unity, in agreement 
with (fl(f). Therefore, the PCPD in the inactive phase confirms the universality of the correlation length finite-size 
scaling amplitude as derived in section 2. Since the entire inactive phase is expected to be in the same universality 
class, that result should apply even to those portions of the inactive phase where our relatively short chains did not 
permit us to carry out a precise extrapolation. 
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2. Second, we briefly discuss the annihilation-coagulation model (]42j). It is well known |£3|.|G0j |G3| that the quan- 
tum Hamiltonian H = H(a,j) is similar to the quantum Hamiltonian H(a + 7,0) of pure annihilation. Therefore 
eqs. ( f45| , fl6| ) also apply to this model, independently of a and 7. If we take a + 7 = 1, the steady-state particle density 
amplitude limi—nx, Lp(L) = (1 + a/j)/(l + 2a/ j) depends on the branching ratio a/7 and is not universal. 

Eqs. (^5j,^6j) state that for the finite-size amplitude of T is independent of the ratio r = a + 7 of the reaction rate 
and the diffusion rate d. In the light of the universality hypothesis of sections 2 and 3, the observed r-independence 
of the amplitude L 2 T means that the critical exponents of the pair annihilation process 2A — > (or the equivalent 
coagulation process 2 A — > A 53 6(J) should also be independent of r, i.e. the mean particle density p(t) ~ t~ 5 with 
5 = 1/2. While that had been anticipated long ago by many people, exact lattice calculations only exist for r = 1, 
see p0| , [47|| . T he only published verifications of the r-independence of 6 we are aware of either used purely numerical 
methods|^!4]|5|], a real -space renormalisation group scheme [ |66| or other ad hoc approximations ]67| . 

The steady-state particle density should scale as 



PL 



D C 2 L-^ Y' (0, C s 



hL 



d+e-0/v± 



h=0 



(50) 



(borrowing a notation from section 2, where the prime indicates the derivative with respect to the second argument and 
h parametrises a source of particles). That is standard finite-size scaling without any readily identifiable universal 
amplitude and is therefore weaker than the forms of sections 2 and 3. It is known from field theory that C2 is 
independent of r |54| while it does depend on the branching ratio a/7 [|33| which is compatible with our results. 

3. Having checked the scaling function universality in some examples with known behaviour, we now illustrate how 
this universality might be used as a diagnostic tool. Appealing to the experience gathered in equilibrium systems [0, 
universal amplitudes might be expected to vary considerably more between distinct universality classes than critical 
exponents. Therefore, even an approximate determination of universal amplitude may allow to conclude on the 
universality class of the model at hand. Reconsider the phase diagram of the pair contact process in fig. |l| Presently, 
there is no consensus on how many universality classes should be realised along the transition line between the active 
and inactive phases for d 7^ 0, 1. As a starting point, we might consider pair mean field theory, which predicts two 
distinct universality classes along the two segments of the pair mean field transition curve (39| . The calculation of 
both steady-state and time-dependent critical exponents in ID from simulations p6[ ] appear to be in agreement with 
this prediction. On the other hand, DMRG studies [[59| and different simulations in ID Q,^,^8| only find evidence 
for a single universality class along the transition line. 

Here, we consider the ratio R = E3 /E2 of the two lowest non- vanishing eigenvalues of H for free boundary conditions 
as obtained from the DMRG. We use the Monte Carlo estimates for p c (d) obtained by Grassberger Ej|. R is equal to 
the ratio of two distinct relaxation times and from eq. (|2~4|), we expect R to be constant within a given universality 
class. Our numerical results are shown in fig. 0. 




FIG. 2. Ratio R — E3/E2 of the two lowest eigenvalues of H at the critical line p = p c {d) of the pair contact process, for 
free boundary conditions and several values of d. 



At first sight, it might be possible to separate the values of R into two classes, one for smaller values of d (up to 
d ~ 0.4 — 0.5) and a limit R^ ~ 2, and one for larger values of d (above d « 0.8) with a limit R^ ~ 3 — 4. The fact 
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that Roo is quite independent of d for 0.2 < d < 0.4 confirms the expected universality (Roc = 1 in the inactive phase 
because there all levels are two-fold degenerate, see appendix). 

However, closer inspection reveals that already for d = 0.5, the data for Rl for L small start out close to 4 and 
then begin to cross over to values close to 2. That signals the presence of strong transient effects in this model. In 
addition, for the values of L for which data exist, the values of R for both d = 0.8 and d = 0.9 are close together and 
quite close to 4. However, these data, in particular for d = 0.8, also show indications that they might also cross over 
to smaller values of R if L could be increased further. The lattice sizes available are not large enough to be able to 
distinguish clearly between the possibilities of a single |39 44 4^j6^] or two Q transitions, although the possibility of 
a single transition appears more likely. If there is a change of the universality class along the critical line at all, fig. ^ 
suggests that it certainly should occur for d > 0.5. 



V. CONCLUSIONS 

In this paper, we have tried to generalize Privman-Fisher Q universality to (steady-state) phase transitions with 
dynamical scaling. For equilibrium systems the standard arguments carry over, the main ingredients being translation 
invariance and hyperscaling. With respect to more standard renormalisation group arguments, only the absence of 
non-universal metric factors in front of the scaling function for £±j in eq. (|23|) is new. Equivalently, this may be 
stated, see eq. ([24T), as the universality of all ratios ,j- 

Out of equilibrium, new arguments must be sought. In the special case of directed percolation, the peculiar 
properties eqs. (|33[ , |35| ) were seen to be sufficient for amplitude universality. These properties may or may not be 
available in other universality classes, but we found some numerical evidence in several reaction-diffusion systems 
that the universal finite-size scaling forms (|23| , p4| ) might indeed hold in general. In the annihilation-coagulation model 
and the pair contact process (inactive phase), we found that the finite-size scaling amplitude of the leading relaxation 
time is indeed independent of the irrelevant parameters we considered. Further evidence in favour of universality was 
found by studying the active-inactive transition line in the pair contact process. In addition, our data for the ratio of 
the two leading relaxation times appear to favour a single universality class along that transition line for < d < 1. 
It remains an open question how to derive the universal scaling forms (|2^ , p4] ) in general non-equilibrium systems. 

Recently, the reaction-diffusion process 2 A — » A and A$A — » 3 A has been studied ]69[ . If the coagulation rate is 
equal to the diffusion rate, the model is exactly solvable. It remains in the universality class of diffusion-annihilation 
for all values of the particle production rate A. For periodic boundary conditions, the exact amplitude of the leading 
inverse relaxation time is limx^oo L 2 Yl — 2ir 2 d, independently of A, as to be expected from ( ^Jitf ) . 

Finally, upon identification the universal finite-size scaling amplitude of the transverse 'spatial' correlation lengths 
£_[_,», it might be tempting to ask if, in analogy to equilibrium P-[To|, there could be relations of universal finite-size 
scaling amplitudes with some exponents. To answer this question would require a set of worked-out examples on 
which some hypothesis of this kind could be tried out. 
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APPENDIX: BETHE ANSATZ CALCULATION 

We calculate the lowest gap T = E2 for p = 1 in the pair contact process and derive the amplitude a in (Q). 

In the p = 1 case, only the pair annihilation AA — > 00 survives. It is well known [j52| that in this case, the quantum 
Hamiltonian can be decomposed H = Hq + Hi in such a way that the eigenvalue spectrum of H is independent of Hi , 
viz. spec (H) = spec (Hq). The latter is related to the spectrum of the XXZ Heisenberg chain which can be found 
from the coordinate Bethe ansatz |5(||37j . 

We begin with the periodic case. The spectrum-generating part of H is 

H a = dH xxz (A,t) + ^(l + d)L (51) 
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where 



1 L 



2 



and 



3d - 1 1 - d 

= 2d~ ' t = d~ ^ 

Since the total spin S z :— J2i=i °f, commutes with Hxxz, the eigenstates Ho\tp n ) = E n \tp n ) can be classified in terms 
of the number n of reversed spins, viz. S z \il) n ) = (L — 2n)\ip n ). The lowest states with n = 0, 1 have zero energy and 
correspond to the two steady states of the model. The lowest gap T is found in the sector n = 2. From |57|] , one has 

£2 = 2rf(2 - cos k - cos fc') , Lk = 2nl-0(k,k') , Lk' = 2nl' - Q(k' ,k) (54) 

where 

^„ ,/s r, Asin((fc - k')/2) 

Q(k, k') = 2 arctan — , - " - — %f — — (55) 

v ; cos((k + k')/2)-Acos((k-k')/2) v ; 

and /, I' = ±f , ■ ■ ■ are distinct half-integers. The total momentum of the corresponding state is P = k + k' = 
2w(I + I')/L' . The lowest energy gap will have P — 0, or k! — —k. Furthermore, the lowest energy state corresponds 
to the choice I = —I' = 1/2, as can be checked by considering the special case A = (because of the symmetry 
between k and k', all levels are two- fold degenerate). We find 

„ ,. , . Lk — n Asinfc 

r = 4d(l-cosfc) , tan = (56) 

v ; ' 2 1 - A cos k y ' 



For L large, the solution of the second equation (]5q) is 



+ (57, 



L \ 1-AL 

Inserting this into (|5^), we arrive indeed at the first case of eqs. ( |45| , |46[ ) . 

Second, we consider free boundary conditions. The spectrum-generating part of H is 



Ho =dffxxz(A,t,r) + i(l + d)(L-l) (58) 



ff xxz (A,<, r) = -\ \j2 [of of +1 + af<jf +1 + Aa?a? +1 ] +r{af + 0%) + t ^ of 1 



(59) 



where we used ( |5^ ) and r = — 1/2. Again, V will be the lowest energy in the n = 2 sector. From the Bethe ansatz 
p6[ , one has 

E 2 = 2d{2 - cos k - cos k') , e 2i ^ k = ( ^^-) ' e -i®(*.*')+ie(-M') ( 60 ) 

V Mr, A) y 

where fk(a,b) = a — b + e lk . A similar equation holds for k' , where fe and fc' are exchanged with respect to (p0|). The 
lowest excitations in the n = 2 sector are found for k' = (up to a two- fold degeneracy). Since r — A = —1, using 
([55|) and taking the logarithm, we find 

/A fc\ 

r = 2d(l - cos A:) , Lfe = tt(I + 1) - 2 arctan f — tan - J (61) 

where I = 0, 1, 2, 3, . . .. The lowest gap is obtained for / = 0, as can be checked for A = 0. In analogy to the periodic 
case, we rewrite the second equation in ( |6l| ) in the form tan((Lfe — 7r)/2) = — A/(l — A)tan(fc/2). For L large, the 
solution is 
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^i( 1 -I^Al + -) < 62 ' 

from which the second case in eqs. (^H[^6|) follows. 

Evidently, for A = 0, that is d — 1/3, we recover for both boundary conditions the well-known results found from 
free-fermion methods. A second closed solution exists for d = 1/2 and free boundary conditions. Then A = 1/2 and 
the second of eqs. ( [5l|) reduces to (L + l)k — n(I + 1). The lowest gap is obtained for 7 = and we arrive at (|47|) . 
We also remark that for the A — 1/2 XXZ chain with boundary terms such that U q (sl(2)) invariance holds, the Bcthe 
ansatz equation can be solved analytically |pq| . 
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TABLE I. Estimates for the normalized amplitude = L 2 T l / '{2tt 2 d) for several values of p and d and periodic boundary 
conditions. The line labeled oo gives the L — > oo extrapolations obtained from the BST extrapolation algorithm J59] and 
uj is the effective correction exponent used in these extrapolations. In some cases the sequences are not monotonous, this is 
indicated by a * at the value of lu used. 







P 


d = 0.10 


L 


0.16 


0.25 


0.50 


0.75 




6 


0.7274475706 


1.1014160997 


1.5543363873 


1.7166144988 




8 


0.6791992625 


0.9951704118 


1.3575537598 


1.4819774944 




10 


0.6559859469 


0.9505431819 


1.2648992305 


1.3631372237 




12 


0.6480035923 


0.9318542215 


1.2110759662 


1.2913069885 




14 


0.6483265043 


0.9247019281 


1.1757710267 


1.2431938517 




16 


0.6530990879 


0.9229800785 


1.1507547386 


1.2087157132 




18 


0.6602551037 


0.9238966184 


1.1320630670 


1.1827963304 




20 


0.6686545657 


0.9260886593 


1.1175466079 


1.1626011550 




oo 


1.04(3) 


1.00(3) 


0.9999(1) 


1.0001(2) 


LJ 


1* 


1.11* 


1 


1.0032 



d = 0.30 


L 


0.20 


0.30 


0.50 


0.75 




6 


0.3393736439 


0.5565119224 


0.8448999662 


1.0205743378 




8 


0.3708396385 


0.6016754416 


0.8664772095 


1.0164775547 




10 


0.4001274702 


0.6385279575 


0.8836745934 


1.0137328840 




12 


0.4287823558 


0.6698979484 


0.8975656827 


1.0118226252 




14 


0.4562391902 


0.6967097472 


0.9088380138 


1.0104117821 




16 


0.4820573931 


0.7197719798 


0.9180744129 


1.0093198954 




18 


0.5060692639 


0.7397682858 


0.9257311017 


1.0084452275 




20 


0.5282757026 


0.7572450703 


0.9321538891 


1.0077262607 




oo 


0.96(4) 


1.0001(1) 


1.0000(3) 


0.99999(3) 


LJ 


1.085 


1.074 


1 


1 



d = 0.50 


L 


0.30 


0.40 


0.50 


0.99 




6 


0.2634463944 


0.3594043077 


0.4362871259 


0.5950441871 




8 


0.3071022508 


0.4156753699 


0.4975263719 


0.6582403419 




10 


0.3480300086 


0.4649432503 


0.5485668186 


0.7054235500 




12 


0.3861353468 


0.5076282318 


0.5909715317 


0.7416603392 




14 


0.4210763888 


0.5445274246 


0.6264865185 


0.7702258735 




16 


0.4528601529 


0.5765616531 


0.6565615085 


0.7932568153 




18 


0.4816890281 


0.6045559890 


0.6823105808 


0.8121848167 




20 


0.5078408510 


0.6291933618 


0.7045790759 


0.8279970372 




oo 


1.00(2) 


1.007(2) 


0.996(2) 


1.0002(2) 


LJ 


0.98 


1 


1.067 


0.994 



d = 0.90 


L 


0.50 


0.75 


0.90 


0.95 




6 


0.0541722289 


0.0729000818 


0.0782390633 


0.0790103464 




8 


0.0671363692 


0.0907802309 


0.0975809620 


0.0985881954 




10 


0.0801954402 


0.1084170765 


0.1165093810 


0.1177191004 




12 


0.0931636953 


0.1256311545 


0.1348893588 


0.1362791328 




14 


0.1059635357 


0.1423692854 


0.1526874700 


0.1542393436 




16 


0.1185589168 


0.1586177799 


0.1699027666 


0.1716014095 




18 


0.1309311714 


0.1743785534 


0.1865468998 


0.1883788379 




20 


0.1430702131 


0.1896605667 


0.2026368456 


0.2045900568 




oo 


1.01(1) 


0.999(2) 


1.001(3) 


1.001(2) 


LJ 


1.081 


1.08 


1.078 


1.079 
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TABLE II. Estimates for the normalized amplitude = L 2 T l / '{tt 2 d) for several values of p and d and free boundary 

conditions, and the L — * oo extrapolation. When the dmrg algorithm did not yield stable results the corresponding finite-size 
entries are left empty. 







d 


p = 0.40 


L 


0.333 


0.50 


0.666 




12 


8572577379 


6725747038 


4583935580 




14 


8737590556 


70538407Q0 


4981 81 5051 




16 


0.8869945040 


0.7322840138 


0.5327541136 




18 


0.8978251485 


0.7547785176 


0.5630353770 




20 


0.9068361253 


0.7738292422 


0.5897524204 




22 


0.9144401079 


0.7901798864 


0.6134826888 




24 


0.9209358777 


0.8043583176 


0.6346890106 




26 


0.9265447049 


0.8167680285 


0.6537455074 




28 


0.9314335748 


0.8277177891 


0.6709574987 




30 


0.9357318117 


0.8374491423 


0.6865761177 




32 


0.9395733083 


0.8461522447 


0.7008078882 




34 


0.9429388328 


0.8539906968 


0.7138314424 




oo 


0.999(3) 


0.9999(5) 


0.9995(9) 




LO 


1.08 


1.0245 


1.029 




p = 0.50 


L 


0.333 


0.50 


0.666 




12 


0.9127350184 


0.7395843960 


0.5254151900 




14 


0.9239698644 


0.7686830752 


0.5656063523 




16 


0.9327349799 


0.7919832846 


0.5996949840 




18 


0.9397547158 


0.8110625710 


0.6289496021 




20 


0.9454946130 


0.8269700162 


0.6543147959 




22 


0.9502701952 


0.8404325820 


0.6765067894 




24 


0.9543020404 


0.8519708198 


0.6960780276 




26 


0.9577542440 


0.8619676543 


0.7134610277 




28 


0.9607290383 


0.8707097768 


0.7289990177 




30 


0.9633226275 


0.8784202612 


0.7429679437 




32 




0.8852717697 


0.7555906209 




34 




0.8913924472 


0.7670514417 




oo 


0.9996(6) 


0.9996(5) 


0.9996(5) 




LO 


1.0925 


1.0165 


1.0191 







p 


d = 0.90 


L 


0.60 


0.70 




0.80 


0.90 




12 


0.1863167451 


0.2042528770 




0.2172237514 


0.2255705980 




14 


0.2119456015 


0.2317569546 




0.2459600505 


0.2549873770 




16 


0.2362199311 


0.2576333440 




0.2728737999 


0.2824706979 




18 


0.2592237747 


0.2820034649 




0.2981179611 


0.3081867029 




20 


0.2810372384 


0.3049794681 




0.3218249399 


0.3322848660 




22 


0.3017367417 


0.3266645014 




0.3441189613 


0.3548989577 




24 


0.3213939696 


0.3471482795 




0.3651117354 


0.3761534335 




26 


0.3400759780 


0.3665336490 




0.3849046790 


0.3961552283 




28 


0.3578455055 


0.3848833645 




0.4035899034 


0.4150046251 




30 


0.3748295723 


0.4022778186 




0.4212504431 


0.4327914966 




32 


0.3908774672 


0.4187863074 




0.4378639789 


0.4495863889 




34 


0.4062441268 


0.4344650723 










oo 


0.986(20) 


0.98(3) 




0.983(15) 


0.984(15) 


to 


1.09 


1.103 




1.100 


1.099 



1G 



